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CHAMBERS  FOR  THE  ELECTRODYNAMIC  CONTAINMENT 
OF  CHARGED  AEROSOL  PARTICLES 


I.  INTRODUCTION. 

Under  proper  conditions  an  electrically-charged  aerosol  particle  can  be  contained  at 
a point  in  space  by  an  alternating  electric  field.  This  phenomenon  may  be  used  to  make  such  a 
particle  accessible  for  a variety  of  physical  measurements.  The  obvious  advantage  of  electrical 
suspension  is  to  avoid  artifacts  required  by  mechanical  suspension. 

A positively-charged  particle  would  remain  at  rest  in  stable  equilibrium  only  at  the 
point  of  minimum  electrical  potential  within  a potential  well.  However,  Eamshaw’s  Theorem 
states  that  there  cannot  be  an  electrostatic  potential  well  in  free  space;  therefore,  containment 
via  electrostatic  fields  of  a charged  particle  at  a point  in  free  space  is  impossible. 

On  the  other  hand,  under  the  proper  conditions  an  alternating  electric  field  permits 
containment  by  producing  on  a charged  particle  a net  time-averaged  force  directed  toward  a 
single  point  in  space.  The  physical  basis  underlying  this  phenomenon  may  be  explained  by 
visualizing  an  electrically-charged  particle  oscillating  in  an  alternating  electric  field.  If  this  field  is 
uniform  (as  between  two  flat,  parallel,  conducting  plates),  the  charged  particle  oscillates  at  the 
driving  frequency  of  the  ac  voltage,  and  the  net  (time-averaged)  force  is  zero.  If,  however,  an 
ac  voltage  is  applied  to  an  electrode  with  an  arbitrarily  curved  geometric  shape,  a nonuniform 
electric  field  results,  and  an  oscillating  charged  particle  experiences  net  (time-averaged)  forces  in 
the  directions  of  decreasing  magnitude  of  instantaneous  field  strength.  A charged  particle  can  be 
contained  within  a small  region  of  space  by  an  alternating  voltage  on  an  electrode  configuration 
whose  electric  field  has  a constant  null  point,  e.g.,  on  a ring  electrode.  The  net  (time-averaged) 
force  toward  the  null  point  is  used  to  counteract  gravitational  or  other  constant  forces. 

H.  Strauble,*’^  using  a single-ring  electrode  connected  to  a high-voltage,  low- 
frequency  (SO-hertz)  source,  succeeded  in  suspending  electrically  charged  aerosol  droplets  with 
diameters  between  30  and  100  pm.  R.  F.  Wuerker  et  al.^  were  able  to  contain  small,  electrically 
charged,  metallic  spheres  in  an  evacuated  chamber  using  hyperbolic  electrodes  and  gave  a brief 
theoretical  discussion  on  the  general  motion  of  the  contained  particles.  A.  Muller^  gave  a more 
detailed  theoretical  treatment.  For  more  recent  applications  of  the  same  principles,  the  reader  is 
referred  to  the  work  of  K.  S.  Fansler  et  al.,^  J.  W.  Schweizer  and  D.  N.  Hanson,*  S.  Ataman  and 
D.  N.  Hanson,''  and  T.  G.  Owe  Berg  and  T.  A.  Gaukler.® 

Our  primary  objective  was  to  determine  the  best  electrode  design  for,  and  the 
usefulness  of,  the  phenomenon  of  electrodynamic  containment  as  an  aid  in  basic  studies  of 
aerosol  particles.  To  achieve  this  objective,  theoretical  calculations  for  idealized  geometries  were 
compared  to  experimental  measurements  for  two  electrode  configurations. 

For  the  purpose  of  computation  it  was  assumed  that  the  electric  field  varied  linearly 
with  distance  in  each  direction  from  the  null  point,  the  vector  proportionality  constant  being  a 
function  of  direction;  this  assumption  resulted  in  a linear  differential  equation  of  motion  readily 
transposable  into  a modified  Mathieu  equation.  The  numerical  calculations  were  performed  using 
standard  methods. 
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Two  electrode  configurations  were  chosen  for  the  experimental  work.  A bihyper- 
boloidal  electrode  configuration  had  the  advantages  of  producing  a radially  linear  electric  field 
strength  while  completely  enclosing  the  suspended  particle.  The  chamber  with  the  bihyper- 
boloidal  configuration  made  possible  accurate  verifications  of  the  theoretical  results.  It  was 
further  equipped  to  operate  under  widely  varying  internal  ambient  conditions.  As  a result,  this 
system  has  become  a sophisticated  tool  for  use  in  aerosol  research.  The  second  configuration 
consisted  of  a ring  electrode  and  two  hemispheric  electrodes.  Although  its  resulting  field  strength 
is  only  approximately  radially  linear,  the  simplicity  of  construction  of  this  system  was 
appropriate  for  initiating  the  work.  Both  configurations  are  described  in  detail  in  the  next 
sections. 

II.  DESCRIPTION  OF  ELECTRODE  CHAMBERS  AND  ASSOCIATED  APPARATUS. 

The  bihemispheric  system,  the  bihyperboloidal  system,  the  associated  circuitry,  and 
the  apparatus  for  inserting  charged  particles  are  described  in  the  following  sections. 

A.  The  Bihemispheric  Electrode  Configuration. 

Figures  1 and  2 show  the  apparatus  for  this  case.  Figure  1 is  a photograph  of  the 
polystyrene  and  wood  chamber.  The  brass  ring  electrode  was  connected  to  the  ac  high-voltage 
power  supply  through  a voltage-dividing  circuit  that  supplied  the  necessary  voltage  to  the  ring 
(49.4%  of  the  ac  voltage  on  the  hemispheres)  to  insure  an  approximately  linear  field  strength 
around  the  null  point.  Removable  microscope  slides,  used  as  illumination  and  observation 
windows,  allowed  introduction  of  the  aerosol  particles.  The  ring  also  had  a small  arc  removed  to 
allow  access  of  the  particle-entry  track  to  the  interior  of  the  chamber. 

Figure  2 shows  the  schematic  diagram  of  the  electrical  circuitry  for  the  bihemi- 
spheric configuration.  The  variable  ac  voltage  unit  could  be  replaced  with  a variable-frequency 
device  for  obtaining  driving  frequencies  other  than  60  hertz.  The  dc  voltage  supply  was  designed 
to  apply  equal,  opposite  polarity  voltages  to  the  two  hemispherical  electrodes,  leaving  the  ring 
plane  at  zero  potential;  the  static  field  could  then  be  calculated  as  that  of  the  idealized  geometry. 
The  continuously  variable  dc  power  supply  was  capable  of  placing  a maximum  potential  difference 
of  200  volts  between  the  hemispheres. 

One  source  of  particle  illumination  was  a strobe  light  positioned  behind  the  chamber 
and  facing  the  observer.  When  the  strobe  light  was  flashing  at  the  same  rate  as  the  ac  driving 
frequency,  a single  sharp  image  of  the  particle  was  visible  through  the  microscope.  With  this  type 
of  illumination  particle-size  measurements  could  be  made  using  a micrometer  FDar  eyepiece  in 
the  microscope. 

An  alternative  light  source  consisted  of  a strip  filament  lamp  with  a vertical  beam 
focused  on  a point  in  the  vicinity  of  the  null  point  from  below  the  axis  of  the  system.  With  this 
device,  the  moving  particle  could  be  viewed  as  a short  white  line  against  a much  darker 
background. 
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Figure  2.  Circuitry  for  Bihemispheric  Configuration 


B.  The  Bihyperboloidal  Electrode  Configuration. 

1 . General. 

Figures  3,  4,  and  5 depict  the  apparatus  for  this  case.  Figure  3 is  a photograph  of 
the  Teflon  chamber,  showing  the  chrome-plated,  brass  electrodes.  The  chamber  was  thermally 
insulated  by  placing  it  entirely  inside  a cube  of  styrofoam.  Figure  4 is  a cross-section  of  the 
chamber. 


To  produce  the  desired  electric  field,  the  electrodes  must  be  surfaces  described  by 

the  equation 

2z2  - = C± 


1 


where  z and  r are  shown  in  figure  4 and  C±  are  constants,  a positive  constant  for  the  two-sheet  j 

hyperboloid  forming  the  top  and  bottom  of  the  chamber,  and  a negative  constant  for  the  j 

one-sheet  hyperboloid  forming  the  sides  of  the  chamber.  j 

These  electrodes  produce  an  electric  field  strength  that  is  zero  at  the  origin  and  has  1 

the  described  linearity,  provided  that  the  two  sheets  are  at  the  same  potential,  as  shown  in  I 

section  HI.  I 
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The  top  and  bottom  hyperboloidal  sheets  were  separated  from  the  central,  single- 
sheet hyperboloid  by  two  horizontal  Teflon  rings  and  could  be  fastened  rigidly  to  the  central 
hyperboloid  with  a set  of  six  plastic  screws  and  nuts.  In  practice,  only  the  bottom  sheet  was 
fastened;  the  top  sheet  simply  rested  in  place  for  easy  removal. 

Two  collinear  holes,  1/8  inch  in  diameter,  were  drilled  diametrically  opposite  each 
other  at  the  narrowest  portion  of  the  electrode.  In  viewing  the  suspended  particles,  one  hole, 
sealed  from  the  outside  by  a cylindrical  glass  window,  formed  an  illumination  port;  the  second 
hole,  a viewing  port,  formed  an  optic  aperture  for  the  microscope  objective  lens. 

2.  Electrical  System. 

The  circuitry  used  to  produce  the  particle-suspending  alternating  electric  field  and 
the  balancing  constant  electric  field  is  shown  in  figure  5.  The  ac  potential  was  impressed  upon 
the  one-sheet  hyperboloidal  conductor.  The  dc  potential  was  impressed  symmetrically  on  the  top 
and  bottom  sheets  of  the  two-sheet  hyperboloidal  electrode,  above  and  below,  respectively,  the 
common  ground  with  respect  to  which  tiie  ac  potential  oscillated. 

The  symmetry  of  the  dc  potential  was  maintained  by  a series  of  variable  potenti- 
ometers for  a range  of  potential  differences  up  to  200  volts.  A scaling  switch  permitted  the  most 
accurate  use  of  the  two  synchronized  helipots  that  controlled  this  voltage. 

The  ac  supply,  variable  up  to  4,000  volts,  originated  from  a high-voltage  trans- 
former. The  input  to  this  transformer  could  be  either  from  line  voltage  or  from  a Belman-Invar 
Invertron  power  amplifier  connected  to  a variable  oscillator.  Because  the  study  utilized  a 
60-hertz  ac  electric  field,  line  voltage  was  acceptable  as  input. 

A Kintel  Model  301  standard  dc  voltage  supply,  a potentiometer,  and  precision 
resistance  boxes  were  used  to  calibrate  both  ac  and  dc  potentials.  The  voltages  were  monitored 
with  a Honeywell  model  1108  Visicorder  and  a Tektronix  type  564  storage  oscilloscope.  The 
dc  voltage  could  be  determined  to  an  accuracy  of  ±0.5  volts,  and  the  ac  voltage  to  an  accuracy  of 
±10  volts. 


3.  Observation  Systems. 

The  two  modes  of  observing  suspended  particles  were  direct  observation  through  a 
microscope  eyepiece  and  photographic  observation  utilizing  a 32-millimeter  movie  camera. 

Lighting  for  both  of  these  modes  was  provided  by  a General  Radio  type  1538-A 
stroboscope  firing  through  the  back-lighting  port.  The  strobe  frequency  could  be  manually  varied 
from  110  to  150,000  flashes  per  minute.  Because  the  apparatus  normally  operated  at  60  hertz, 
the  stroboscope  was  automatically  synchronized  with  the  60-hertz  line  voltage.  A time-delay 
device  permitted  the  observation  of  particles  at  any  point  in  their  60-hertz  vertical  oscillations. 
A circuit  that  automatically  doubled  the  stroboscope  frequency  to  1 20  hertz  was  also  introduced 
to  produce  double  images  of  the  particles.  With  the  proper  time  delay  the  separation  of  the  two 
images  of  a particle  would  reach  a maximum  corresponding  to  the  oscillation  endpoints 
separation. 


The  particle  image  seen  through  the  second  port  was  formed  using  a Bausch  and 
Lomb  32-mm  focal  length  microscope  objective  lens  with  a numerical  aperture  of  0.10  in 
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coixjunction  with  a Vickere-AEI  image^plitting  eyepiece  mounted  at  the  end  of  a 6-inch -long 
optical  tube  equipped  with  an  adjustable  iris.  The  optical  apparatus,  supported  by  a tripod,  was 
physically  independent  of  the  chamber.  When  viewing  the  particle,,  the  observer  adjusted  the 
eyepiece  tube  to  touch  the  tube  containing  the  objective  lens,  obtaining  a focused  image  with  a 
fixed  magnification.  The  eyepiece  could  then  be  calibrated  for  measuring  particle  diameters.  Its 
image-splitter  consisted  of  a prism  controlled  by  a ruled  micrometer.  By  bringing  the  two  split 
images  into  apparent  contact,  a diameter  measurement  could  be  made  in  the  range  10  to  250  ;xm 
with  an  accuracy  of  ± 1 pm. 

The  32-mm  movie  camera  required  a bright  source  of  light  flashing  at  a suitable 
frequency  while  the  camera  shutter  was  open.  These  requirements  were  met  by  use  of  a multiple 
timer,  additional  electronic  circuitry,  and  a large  capacitor  attached  to  the  stroboscope.  The 
resulting  film  records  of  suspended  particles  were  useful  in  analyzing  the  changes  in  shape  and 
size  of  the  particles. 

4.  Thermal  System. 

The  top  and  bottom  sheets  of  the  bihyperboloid  electrode  configuration  contained 
reservoirs  for  liquid  (see  figure  4)  that  were  used  to  control  the  temperature  of  the  chamber. 
Once  a given  chamber  temperature  had  been  established,  it  was  maintained  by  the  insulating 
foam  enclosure.  Because  measurements  of  the  dc  balancing  voltage  and  the  ac  splitting  voltage 
(see  section  III)  were  severely  affected  by  thermally-induced  convection  currents,  the  observer 
had  to  wait  for  thermal  equilibrium  for  1/2  hour  from  the  time  liquid  circulation  ceased  before 
taking  measurements. 

A small  thermocouple,  inserted  through  the  lower  Teflon  ring  and  connected  to  the 
Honeywell  Visicorder  allowed  the  temperature  of  the  interior  of  the  chamber  to  be  determined 
to  within  ±rc. 

5.  Atmospheric  System. 

Ports  in  both  the  top  and  bottom  Teflon  rings  were  used  to  control  the  ambient 
atmosphere  in  the  chamber.  With  appropriate  equipment  connected  to  these  ports,  the  desired 
pressure  and  composition  of  the  gas  in  the  chamber  could  be  established. 

In  the  calibration  experiments,  very  low  ambient  humidity  was  obtained  by  slowly 
circulating  air  first  through  phosphorus  pentoxide  and  then  through  the  chamber.  This  air 
circulation  was  stopped  during  measurements. 

6.  Particle  Insertion  Device  and  Method. 

Figure  6 is  a photograph  of  the  apparatus  used  to  insert  aerosol-sized  particles  into  a 
chamber.  This  apparatus  consisted  of  three  22-inch,  parallel,  metal  rods.  One  rod  was  placed 
below  the  other  two  rods  so  that  the  cross-section  formed  a V-shaped  track  supported  by  a 
small,  notched,  plastic  rectangle  on  one  end  and  by  a larger  rectangular  plastic  shield  on  the 
other  end.  The  top  two  rods,  joined  behind  the  plastic  shield  by  a loop,  were  connected  to  the 
ac  voltage  source.  The  bottom  rod  was  rounded.  The  open  end  of  the  track  was  placed  in  the 
chamber  (in  the  case  of  the  bihyperboloidal  chamber,  after  removal  of  the  upper  sheet). 
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Figure  6.  Aerosol  Particle  Insertion  Apparatus 


The  liquid  or  dissolved  solid  to  be  studied  was  placed  inside  a glass  capillary  tube 
that  was  connected  to  a variable  dc  high-voltage  source  normally  operating  between  4,000  and 
10,000  volts.  A polydisperse  spray  of  charged  particles  was  then  obtained  by  striking  the 
capillary  tube,  causing  it  to  vibrate. 

These  particles  were  allowed  to  fall  through  the  loop  connecting  the  two  upper  rods, 
where  particles  with  appropriate  charge  and  mass,  depending  on  the  voltage  on  the  rods,  were 
captured;  the  track  was  then  tilted  downwards  toward  the  chamber  to  move  the  particles.  All 
particles  except  the  one  selected  for  observation  were  removed  from  the  track  (or,  if  necessary, 
from  the  chamber)  by  means  of  plastic  rods  that  attracted  the  unwanted  particles  either  by 
induction  or  by  virtue  of  being  charged  oppositely  to  the  charge  on  the  particles.  Once  the 
selected  particle  was  in  the  chamber,  the  track  was  removed.  In  the  case  of  the  bihyperboloidal 
system,  the  upper  sheet  was  then  replaced. 

It  was  found  necessary  to  assist  the  initial  capture  of  the  particle  in  the  chamber  by 
impressing  a dc  voltage  on  the  electrodes.  In  the  case  of  the  bihemispheric  system,  this  was 
approximately  a gravity-balancing  voltage.  In  the  case  of  the  bihyperboloidal  system,  however, 
the  ac  field,  distorted  by  the  removal  of  the  upper  sheet,  repelled  the  particle  from  the  apex  of 
the  lower  sheet,  permitting  it  to  slide  down  and  out  of  the  containing  volume;  it  was  therefore 
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necessary  to  impress  upon  the  lower  sheet  a dc  voltage  of  opposite  sign  to  the  gravity-balancing 
voltage  to  attract  the  particle  back  toward  the  apex  of  the  sheet.  This  voltage  was  removed  when 
the  chamber  was  closed. 

IIL  MOTIONS  AND  OBSERVATIONS  OF  CHARGED  AEROSOL  PARTICLES. 


Theoretical  considerations,  discussed  in  the  next  section,  show  that  some  electro- 
dynamic parameters  of  a particle  can  be  used  to  determine  other  physical  properties  of  the 
particle  such  as  charge  and  mass.  These  considerations  are  presented  in  detail  for  the 
bihyperboloidal  system,  for  which  the  equations  of  motion  are  linear  and  hence  reasonably 
tractable.  The  results  for  the  bihemispheric  system  are  essentially  the  same,  differing  only  in 
minor  variations  introduced  by  the  nonlinearity  of  the  ac  field  and,  of  course,  different  values  of 
the  geometric  parameters. 

A.  Solutions  of  the  Equation  r Motion. 

The  equation  governing  the  motion  of  the  particle  in  the  chamber  is 


where 


mR  = -Kj)R  + q(Ej  + E cos  + mg 


R = position  of  the  particle  relative  to  null  point 
m = mass  of  the  particle 

q = charge  on  the  particle  (considered  positive  here) 

Kp  = coefficient  describing  aerodynamic  drag  force 
E cos  wt  = the  alternating  part  of  null-point  field 
Ej  = constant  part  of  null-point  field 
<jj  = angular  frequency  of  applied  voltage 
t = independent  variable,  time 

Edc  = static  field  between  the  upper  and  lower  portions  of  the  outer 
electrode  where  these  parts  are  at  different  potentials 

g = acceleration  due  to  gravity 

NOTE:  Dott  uidicate  dUTeientittion  with  respect  to  time. 

For  small  values  of  the  Reynolds  number.  Re  (=  pd  iRl/ij),  Kp  is  given  by 


Kp  = SiTTjd 


(Stokes’  Law) 
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where  d is  the  particle  diameter,  and  17  and  p are  the  viscosity  and  density  of  the  air,  respectively. 
For  larger  values  of  the  Reynolds  number  a correction  term  must  be  added,  as  in  the  formula  due 
to  Oseen: 

Kj)  = 3irT?d(l  +^Re) 


the  form 


Any  cylindrically  symmetric  potential  function  in  charge-free  space  can  be  written  in 


^ Pn  (cos  ®s) 


where  and  Qg  are  the  radial  and  polar  angle  coordinates  of  a spherical  coordinate  system, 

Pfi  (cos  0^)  are  Legendre  polynomials,  and  the  constant  coefficients  An  depend  on  the  geometry 
of  the  field  source. 

If  the  potential  is  also  symmetric  about  the  plane  Q as  in  the  null-point  fields 

considered  here,  the  odd-term  coefficients  are  all  zero.  If  the  series  is  now  rewritten  in  cylindrical 
coordinates  r,  B,  z the  field  becomes 

^ = Aq  +A2ir'^  - 2z2)  + A/^{3r^  - 24r2  + ^4)  (2) 

+ A(,  (5r6  - 90A2  + I20r2z4  _ I6z6)  + . . . 

where  A ^ are  a new  set  of  constant  coefficients.  In  the  case  of  the  bihemispheric  configuration, 
the  potentials  relative  to  a distant  ground  on  the  ring  and  on  the  spheres  were  adjusted  so  that 
A4  = 0,  leaving  in  the  series  (besides  the  quadratic  term)  only  terms  of  sixth  and  higher  powers, 
which  are  very  small  in  the  region  around  the  origin.  The  potential  is  then  approximated  by 

* = Aq  + A2(/’2  - 2z2)  (3) 

If  the  quartic  and  higher  terms  of  equation  2 are  dropped  altogether,  so  that  equation  3 is  exact, 
the  equipotentials  are  hyperboloids  of  revolution  described  by 


,2.2,2  = *^ 


the  right-hand  side  being  constant  on  an  equipotential.  The  electrodes  in  the  bihyperboloidal 
configuration  are  described  by 

r2  - 2r2  = ro2  (5) 

r2  - 222  = -2zq2  (6) 
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where  tq  is  the  radius  of  the  one-sheet  hyperboloid  and  22q  is  the  distance  between  the  upper  and 
lower  sheets  of  the  two-sheet  hyperboloid.  The  values  of  >4o  and  A 2 can  be  found  by  noting  that 
the  two-sheet  hyperboloid  is  grounded  (♦  ■ 0),  and  that  the  potential  ♦ on  the  one-sheet  hyper- 
boloid is  -K  cos  ojt  (where  K is  the  peak  voltage),  equating  the  right-hand  sides  of  equations  4 
and  5,  and  equating  the  right-hand  sides  of  equations  4 and  6.  When  these  values  are  substituted 
into  equation  3,  the  expression  for  becomes 


In  this  investigation  higher  power  terms  are  dropped  and  is  approximated  as 


Edc  - A'^k-  CAnVAfk 


dc'^dc* 


where  is  the  potential  difference  between  the  sheets  (Fjc  >0  if  the  lower  sheet  is  at  the 
higher  potential)  and  Cjc  is  an  experimentally  determined  geometric  constant.  Equation  9 is 
exact  in  two  cases,  to  wit:  (1)  when  FJjc  = 0 and  the  particle  oscillates  below  the  null  point 
only  under  the  influence  of  the  null-point  ac  field  and  gravity,  and  (2)  when  the  particle  is 
balanced  at  the  null  point  because 

r and  z are  then  zero  and  the  higher  power  terms  of  equation  8 vanish. 

The  principal  situation  of  interest  in  this  investigation  for  which  equation  9 is  not 
exact  occurs  when  the  particle  is  balanced  at  the  null  point  but  is  oscillating  about  the  null  point 
because  it  is  in  an  unstable  region  of  the  k-E  chart  (described  later  in  this  section). 

The  axis  of  symmetry  of  the  apparatus  is  vertical  so  that 

g * -gk 

To  transform  equation  1 into  dimensionless  form,  let 

T = cur  (lOa) 

^1  = (10b) 

^2  = yid  (10c) 

f = z/d  (lOd) 

where  x,  y,  and  z are  the  cartesian  coordinates  of  R,  and  d is  the  particle  diameter.  Let  equa- 
tion 10  be  substituted  into  equation  1 and  the  latter  be  divided  by  Then  equation  1 written 
in  terms  of  its  components  is 

f"  + (k  + Kju)f'  + (£  cos  T + £j)f  + h - 0 (11) 

r + (K +Ki«)f'-^£  cos  r + £i){  = 0 (12) 
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where  the  primes  denote  differentiation  with  respect  to  r,  and  where 


18n 


u » iR'I/d 
^1- 

urm 


h = -L- 


Cdcft 


(ji^d  co^dm 


»"dc 


i represents  (-j  or  ^2"  Pj  is  the  density  of  the  particle  material,  and  C is  the  geometric  constant 
relating  the  potenti^  on  the  electrodes  to  the  null-point  field  strength.  From  equation  7 it  is 
easy  to  see  that  for  hyperboloidal  electrodes, 


C = 


+ 2zo^ 


Because  tq  = 1.25  cm  and  zq  = 1.75  cm  for  the  device  described  in  this  report,  C = 5,203  m" 


If  the  aerodynamic  drag  is  described  by  Stokes’  Law, 


K j = 0 


If  the  drag  is  described  by  Oseen’s  Law, 


Now  equations  1 1 and  12  are  of  essentially  the  same  form;  for  reasons  that  will 
become  clear  later,  equation  1 1 has  more  physical  significance  than  equation  12,  and  the  discussion 
will  concentrate  on  equation  1 1 but  is  applicable  to  both.  Solutions  are  first  discussed  for  the 
case  of  Stokes’  Law,  and  the  effects  of  deviations  from  Stokes’  Law  are  considered  later.  For 
Stokes’  Law  equation  11  is  linear,  viz.: 


r + icf'  + (£  cos  T + £i)f  + h = 0 


(13) 


The  solution  of  equation  13  has  the  form 


^ + Bi2(T)  + P(r) 


(14) 


where  f j and  {'2  independent  solutions  of  the  homogeneous  equation 

S"  + Kf'  + (£  cos  T + £j  )f  = 0 
and 

oe 

p(r). 

« = -«» 


(15) 


(16) 


The  method  of  calculating  the  coefficients  is  shown  in  appendix  A.  The  coefficients  a„  are 
proportional  to  h and  vanish  if  A = 0.  By  means  of  the  transformation 

T=2r  (17a) 


= (17b) 

equation  1 5 may  be  written  as  the  Mathieu  equation 

T"  + (a  - 2q  cos  2r)f  = 0 

where  a and  q are  constants.  Properties  of  f j and  f 2 1*®  discussed  in  terms  of  the 

well-known  properties  of  solutions  to  the  Mathieu  equation.  From  this  comparison  (see 


appendix  A)  it  is  seen  that  the  solutions  f j and  f2  should  be  of  the  form 

tl  - e'  2" 

where  P(t)  is  a bounded  periodic  function  with  period  2ir,  and  7 can  take  one  of  the  forms 

7=0  (19a) 

7 = a + y/  (19b) 

7 = 01  (19c) 
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where  a and  /3  are  real.  It  is  evident  that  a may  be  considered  positive  without  loss  of  generality. 
The  coefficients  A and  B of  equation  14  can  be  calculated  from  the  conditions  at  any  particular 
value  of  r;  in  general  neither  vanishes.  Therefore  f is  stable  (bounded  for  t > 0,  where  t = 0 
represented  the  initial  situation)  only  if  and  {'2  stable,  i.e., 


/?(T)-y»c<0 


whence 


-i?(7)-yK<0 


where  Ri'y)  is  the  real  part  of  7.  This  condition  is  obviously  satisfied  for  the  case 

y = Pi 

and  requires  that  a < ic/2  in  equation  19.  A method  for  calculating  7 is  shown  in  appendix  A. 

A chart  showing  the  nature  of  7 as  a function  of  k and  E,  for  £ j = 0,  is  shown  in 
figure  7.  In  the  space  defined  by  k and  E there  is  a series  of  regions  where  a ^ 0 that  are 
separated  by  regions  where  7 = Pi.  In  the  former  regions  7 is  alternately  given  by  7 = a and 
7 = a + -Lj.  It  is  evident  that  in  regions  where  7 = a the  periodicity  of  f is  2jr,  and  that  in 
regions  where  7 = a + y i,  the  periodicity  of  f is  4ir;  in  regions  where  7 = pi  the  function  is  in 
general  not  periodic.  On  the  boundaries  between  7 = Pi  and  y = a,  a and  p both  vanish,  so  that 
7 = 0.  On  the  boundaries  between  7 = pi  and  7 = 0+  i regions,  7 = y /.  In  the  interiors  of 
the  regions  where  a # 0 are  regions  where  a > these  interior  regions  the  solution  f of 

equation  13  is  not  stable.  The  curves  on  which  a = -L(c  are  the  boundaries  between  regions  of 
stable  and  unstable  solutions  of  equation  1 3 and  are  designated  as  instability  boundaries  or 
quasi-instability  boundaries  in  this  report.  Figure  7 shows  the  first  such  unstable  region  and  part 
of  a second. 

In  region  /?q  in  figure  7,  7 is  real,  but  a never  exceeds  the  function  f is  stable 
everywhere  in  Rq  and  there  is  no  instability  boundary.  The  periodicity  of  f in  Rq  is  2v.  In  the 
region  I\,  y = Pi,  and  f is  stable  and  nonperiodic.  In  regions  R j ^ and  R ^ fj,  y = a + y /,  and 
the  period  of  f is  47r.  In  £ j a < y »c,  and  f is  stable;  in  /?  j o > -yic,  and  f is  unstable. 

The  boundary  between  /?j^  and  /?j  [/  is  the  first  instability  boundary.  In  the  regions  ^2^ 

/?2  {/.  7 = 0,  and  the  period  of  f is  2ff.  The  second  instability  boundary  separates  RjjS'  where 
a > and  f is  unstable.  Between  Rijs  and  /?2  5 lies  the  narrow  region  I2,  where  7 = Pi,  and 
f is  stable  and  nonperiodic.  The  quantity  v,  shown  in  figure  7 are  curves  of  constant  v,  is  given 
by 

V = a --i»c 
2 
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Figure  7.  Instability  Boundaries 
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and  represents  the  rate  of  expansion  (or  damping  for  < 0)  of  f with  increasing  t.  At  the 
instability  boundary,  v = 0. 


The  dependence  of  the  instability  boundaries  on  the  value  of  £ j is  illustrated  in 
figure  8.  The  most  radical  difference  between  the  £|  = 0 curve  and  the  £|  # 0 curves  is  the 
appearance  of  a new  boundary  in  the  lowest  y = a region  for  the  case  of  £ j < 0.  On  the 
horizontal  axis  of  the  chart  £ = 0,  and  equation  1 5 is 

r + f'  -t  £if  = 0 

Then  f is  given  by 

. _ . -jKt  yr  -4rKr  -yt 
S = Ae  e +Be  2 e 

where 

If  £l  > -^K^,  y = 0i;if  0<  El  < y = a<  ^k;  in  either  case  f is  stable.  If  £i  < 0, 

y - OL>  and  f is  unstable.  The  new  stability  boundary  is  the  upper  boundary  of  this 
unstable  region. 

In  the  stable  regions  of  the  k-E  chart  the  functions  f j and  {"2  become  vanishingly 
small  as  t increases;  except  very  near  the  stability  boundary,  f|,  the  more  slowly  decreasing  of 
the  two  functions,  decreases  to  10“^  of  its  initial  value  (or  of  the  initial  value  of  f)  in  a few 
tens  of  cycles.  Equation  14  describing  f can  then  be  represented  by 

f = P(t)  (20) 

Physically  this  means  that  in  the  stable  region  the  particle  motion  is  described  by  p(r)  alone.  It 
has  been  noted  above  that  if  h = 0,  p(t)  = 0.  Hence  if  the  value  of  E^^  is  set  at  mglq,  so  that 
h = 0,  the  particle  will  be  at  rest  at  the  null  point  if  x and  £ lie  in  the  stable  region. 

If  there  is  no  potential  difference  between  the  two  sheets  of  the  hyperboloid,  the 
force  holding  the  particle  against  gravity  is  supplied  by  the  phase  difference  between  the  applied 
field  and  the  particle  motion.  This  can  be  shown  by  substituting  equation  1 6 into  equation  1 3 
and  integrating  with  respect  to  t over  one  cycle  of  cos  t.  The  result  is 


2iraQ£j  + ir(aj  + a_j  )£  + 2irh  = 0 (21) 


Let 


Note  that  = a*^ . 


when  Rf^  and  /„  are  real. 
Then 


P(r)  = ao 


cos  (MT 


(22) 


where  C„  = 2la„l  and  tan  <(>„  = as  shown  in  figure  9.  The  phase  angle  is  the  angle  by 

which  the  nth  component  of  the  particle  motion  lags  the  driving  voltage. 


Figure  9.  Phase  Angle  in  Particle  Oscillation 


If  equation  21  is  written  in  terms  of  equation  22  and  the  physical  quantities,  it 

becomes 

(“o*'!  + *1 ) ■'  * “I:  ° “ 

o)^m  ^ d 

whence  the  required  restoring  force,  mg,  is  given  by 

mg  = -Cq  (KiAq  + 

We  shall  consider  henceforth  only  the  case  £ j = K|  = 

Then 

mg  = — cos  0])^  (23) 


^KCi  cos  01  )d 
0. 
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In  the  case  now  under  consideration,  where  the  constant  term  in  equation  1 1 is  due 
to  gravity  alone,  oq  < 0,  and  the  particle  path  lies  entirely  below  the  null  point.  The  particle 
position  is  approximately 

z = p(T)d  = [uq  + Cl  cos  (cjt  -4>i)]d  (24) 

and  the  field  strength  is  approximately 

E a (CIuqIK  cos  T)k  (25) 

Equation  25  reflects  only  the  oscillation  of  direction  of  the  field  E.  In  fact,  because 
the  field  becomes  weaker  in  the  direction  of  the  null  point  (upward),  the  particle  is  in  a stronger 
field  when  it  is  below  the  null  point  a distance  IflQc/l  than  when  it  is  above.  Equation  23  shows 
that  cos  must  be  negative.  If  cos  = -1,  so  that  0j  = tt,  equations  24  and  25  show  that 
below  the  null  point  the  particle  is  in  the  stronger  part  of  the  field  when  the  field  is  directed 
upward,  and  in  the  weaker  part  of  the  field  when  the  field  is  directed  downward;  the  net  electric 
force  is  upward,  opposite  the  force  of  gravity. 

Figure  10  shows  the  particle  depression  (the  mean  distance  below  the  null  point  of 
the  vibrating  particle)  when  it  is  being  held  up  only  by  the  null-point  ac  field.  The  constant  c is 
defined  by  the  equation 

V ' ^ 

The  quantity  laQcfl  is  the  particle  depression;  for  a frequency  of  60  hertz  it  is  given  by  the  follow- 
ing expression  in  millimeters: 

\aQd\  = 0.06895c 

The  dotted  lines  in  figure  10  are  the  instability  boundaries.  When  k and  E are  such  that  the 
general  solution  described  by  equation  18  is  stable  and  = 0,  maximum  stability  corresponds 
to  the  minimum  particle  depression.  The  electrodynamic  quantities  that  may  be  varied  are  the 
peak  voltage  V and  the  angular  frequency  co.  The  curve  labeled  "E  minimum”  in  figure  10  is  the 
set  of  (k,£)  pairs  for  which  the  particle  depression  is  minimal  when  k is  constant  and  E is  varied, 
i.e.,  when  the  frequency  is  constant  and  the  voltage  is  varied.  If  only  the  frequency  is  varied,  the 
values  of  k and  E for  which  the  particle  depression  is  minimal  fall  on  the  curve  marked  “frequency 
minimum”.  These  curves  are  of  significance  when  the  particle  is  first  being  captured  in  the  chamber 
before  gravity  can  be  overcome  by  applying  the  proper  voltage  between  the  two  hyperboloidal 
sheets;  when  the  proper  voltage  has  been  applied,  the  part  of  the  oscillation  described  by 
p(r)  vanishes  altogether.  The  particle  depression  minimum  curves  can  be  used  to  establish  chamber 
design  requirements  for  containment  of  particles  in  various  size  ranges. 
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Figure  1 


We  have  seen  that  when  Pjc  = 0 the  net  electrodynamic  force  is  proportional  to  the 
cosine  of  0j,  the  angle  by  which  the  fundamental  term  of  the  Fourier  series  p(j)  describing  the 
motion  of  the  particle  after  f j and  f2  are  damped  out  lags  the  driving  voltage  -V  cos  t.  For  small 
values  of  E 


tan  0 3£  -K 


(26) 


A quantity  related  to  the  lag  angle  0j  is  the  value  for  which 

PCtq)  = P(to  + ff) 


(27) 


Experimentally  this  quantity  is  determined  by  observing  the  vibrating  particle  in  the  light  of  a 
stroboscope  flashed  twice  per  cycle  of  the  driving  voltage.  Two  separate  images  of  the  particle 
are  observed  except  when  the  strobe  flash  occurs  at  a point  in  the  cycle  for  which  equation  27 
is  true,  when  the  two  images  coincide.  From  equations  22  and  27 


flQ  cos  (nr  - 0«)  ~ aq  " ^n^ 


«=1 


n=\ 


whence 


^2«+l  [(2n+l)r  - 02;,+i]  = 0 (28) 

n=0 


If  E is  small  enough  the  particle  motion  is  simple  harmonic  and  p(t)  can  be  represented  by 

P(t)  = aq  + Cj  cos  (t  - 0j  ) 

In  this  case  equation  28  becomes 


Cj  cos  (t  - 0i)  = 0 


and  because  if  Cj  vanished  there  would  be  no  electrodynamic  force  to  hold  up  the  particle, 

TO  2 0,  ^ ^ 

Then  from  equation  26, 


cot  tq  a -tan  02  s ic 


27 


r 


For  larger  values  of  E, 


^0 


= -1  + 


C3  C5 

+ — sin  ^3  + 3 — r cos  ^3  sin  'Ps  -jr  s*"  Ps  ^ 


where 


^3  = 3^1  - ^3 


and 


^5  - 5 - 05 


Figure  1 1 shows  0j  and  cot  tq  as  functions  of  k for  several  values  of  E.  Experi- 
mental values  of  cot  tq,  taken  from  observations  on  an  evaporating  particle,  are  also  shown.  The 
value  of  E was  about  0.26  at  the  beginning  of  the  observation  period  when  k was  small,  but 
in  eased  markedly  toward  the  end,  accounting  for  the  large  difference  between  k and  cot  tq 
' ' en  K was  large. 

When  Oseen’s  Law  is  used  to  describe  the  aerodynamic  drag,  equation  11  is  no 
longer  linear,  and  tlie  method  of  continued  fractions  described  in  the  appendix  A is  no  longer 
applicable.  For  this  case  equation  1 1 is  integrated  numerically  using  Taylor  expansions  for  f and 
f'  and  a computer.  The  principal  effect  of  the  velocity  squared  term  is  that,  in  the  unstable 

regions  y and  /?2  u ^'^ure  7,  the  amplitude  of  the  particle  vibration  does  not  expand 

indefinitely  as  required  by  the  solutions  to  the  linear  equation  but  approaches  a constant  ampli- 
tude. Except  near  the  stability  boundary  this  constant  amplitude  is  reached  in  a few  tens  of 
cycles. 


Figure  12  shows  the  amplitude  of  vibration  (the  total  length  of  the  path  of  the 
particle)  as  a function  of  E for  various  values  of  k,  for  the  case  h - 0.  The  curves  are 
labelled  according  to  the  value  of  « to  which  they  correspond.  The  marks  labeled  L0_;(  and 
LI0_j  represent  the  intersections  of  the  x = x line  on  figure  7 with  the  lower  and  upper  boundaries 
of  the  (odd)  instability  region  /?)  y.  LE^  and  UE^^  designate  similar  boundaries  for 
curve  labelled  ic  = 0. 1 represents  a 9’2-pm  particle  at  60  hertz  in  air.  The  amplitude  increases 
from  zero  starting  at  £ = 0.46,  which  is  seen  from  figure  7 to  represent  the  lowest  stability 
boundary.  The  vibration  has  a period  of  47r  (twice  the  period  of  the  driving  voltage)  as  required 
by  the  linear  solution.  The  curve  rises  sharply  from  the  point  E = 0.46  and  continues  past  the 
second  stable  region  (£  = 3.76)  and  into  the  region  R2  increasing  monotonically.  The  perio- 
dicity in  /?2  If  is  4ff,  in  contrast  to  the  periodicity  of  iTr’ required  by  the  linear  solution  in  this 
region.  Since  the  height  of  the  bihyperboloidal  chamber  for  this  case  is  only  190  particle  diame- 
ters, observation  of  this  mode  of  vibration  for  this  particle  size  is  extremely  difficult. 

The  curves  labeled  k = 1 .0  represent  a 29-pm  particle  in  air  at  60  hertz.  The  larger 
of  these  curves  increases  from  zero  at  £ as  1.18  as  required  by  the  linear  solution,  and  the  period 
of  vibration  is  4ir.  Notice  that  this  curve  extends  past  £ = 4.2,  the  upper  boundary  of  R\  ij,  and 
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AMPLITUDE  (PARTICLE  DIAMETERS) 
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Figure  12.  Amplitude  of  Vibration  in  Unstable  Region 
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and  E = 4.4,  the  beginning  o(  R2  (j-  Some  of  the  calculated  points  fell  on  this  extended  portion 
of  this  curve.  Other  calculations  pve  zero  amplitude  for  values  of  E lying  above  the  upper 
boundary  of  j [j  or  lay  on  the  smaller  ic  = 1 .0  curve  rising  from  zero  at  £ = 4.4,  the  beginning 
of  R 2 If.  The  periodicity  of  points  lying  on  the  lower  curve  is  27r,  as  required  by  the  linear 
solution.  Whether  the  amplitude  calculated  for  a given  value  of  £ > 4.2  lies  on  the  larger  or  the 
smaller  (<=1.0  curve  appears  to  depend  at  least  in  part  on  the  initial  values  of  f and  The 
chamber  height  for  this  size  particle  is  600  diameters. 

The  curve  for  k = 2.0,  representing  a 2\-nm  particle,  rises  comparatively  slowly  from 
zero,  the  value  predicted  by  the  linear  solution  it  E = 3.0,  and  increases  monotonically  until 
E = 5.50,  a value  somewhat  greater  than  the  value  E = 5.10  representing  the  upper  boundary  of 
/?j  (/.as  determined  from  the  linear  equation.  The  amplitude  for  E = 5.55,  the  next  higher 
value  for  which  a calculation  was  made,  is  zero.  At  £ = 6.2,  the  lower  boundary  of  £2,!/’ 
amplitude  begins  to  increase  slowly  with  increasing  £;  the  period  of  vibration  is  2jr  on  this  part 
of  the  curve,  as  required  by  the  linear  solutions.  Particles  of  this  size  (20-/im)  and  smaller  can 
readily  be  experimentally  observed  vibrating  in  the  4]r-period  mode  described  here. 

The  instability  curves  for  horizontal  vibration  can  be  seen  from  equation  12  to  be 
the  same  as  for  the  vertical  vibrations  discussed  hitherto,  with  the  following  exceptions: 

( 1 ) For  a given  value  of  V,  the  signs  before  £ and  £ j are  different  from  the  signs 
for  the  vertical  case.  The  sign  of  £ is  immaterial  physically,  since  the  phase  of  t is  arbitrary.  The 
difference  in  the  sign  before  £j  means  that  the  new  instability  region  that  appears  for  vertical 
motion  when  £j  is  negative  appears  for  horizontal  motion  when  £(  is  positive.  Therefore,  for 
any  given  nonzero  value  of  £j,  there  is  an  unstable  region  at  the  bottom  of  the  k-£  chart  either 
for  vertical  motion  or  for  horizontal  motion. 

(2)  The  coefficients  £ and  £ j are  multiplied  by  one-half.  This  means  that 
corresponding  points  on  the  chart  require  twice  as  much  alternating  voltage  for  the  horizontal 
component  of  motion  as  for  the  vertical  component.  Hence  for  values  of  k less  than  about  1.8, 
the  horizontal  and  the  vertical  instability  regions  overlap;  while  for  values  of  k greater  than  1.8, 
there  is  a region  of  horizontal  double-period  oscillation  lying  somewhat  above  the  first  vertical 
instability  region  and  overlapping  the  second  (single-period  oscillation)  vertical  instability  region. 

(3)  With  the  axis  of  the  chamber  vertical,  there  is  no  term  corresponding  to  h. 

Hence  there  is  no  horizontal  equivalent  to  the  part  of  the  motion  described  by  p(t). 

A number  of  phenomena  that  can  be  correlated  with  the  theoretical  discussion  have 
been  observed  experimentally.  Generally,  when  x and  £ lie  outside  the  unstable  regions,  and 
Kdc  “ 0,  the  particle  is  observed  to  oscillate  below  the  null  point  with  the  frequency  of  the 
driving  voltage  (60  hertz  in  this  case),  as  described  by  equation  20.  When  the  weight  of  the 
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particle  is  balanced  electrostatically  by  setting  Fjg  at  the  proper  value,  the  particle  is  observed 
to  come  to  rest  at  the  null  point;  this  corresponds  to  the  fact  that  h vanishes  in  this  case,  and 
because  all  coefficients  of  p(t)  are  proportional  to  h,  p(t)  also  vanishes. 

The  charge-to-mass  ratio  of  the  particle  is  found  by  electrostatically  balancing  the 
particle  at  the  null  point;  it  is  given  by 


± = g 
^ ^dc^dc 

wliere  is  the  voltage  necessary  to  balance  the  particle.  The  value  of  E for  any  particular 
setting  of  the  alternating  voltage  V is  then 


E = 


V 


dc 


The  value  of  has  not  been  calculated  geometrically  for  the  hyperboloidal  chamber;  it  was 
found  experimentally  to  be  17.7  ± 0.1  meter”*,  as  described  later  in  this  report.  Then  for  any 
particle  for  which  Fjjc  has  been  measured,  E corresponding  to  the  alternating  voltage  is 

E = 0.0210—^ 

^dc 

When  particles  of  any  diameter  over  1 9 pm  were  balanced  at  the  null  point  and 
made  to  oscillate  by  increasing  the  alternating  voltage  V,  the  oscillation  was  invariably  in  the 
double-period  mode  required  by  the  theory,  and  the  values  of  experimental  quantities  were 
consistent  with  the  values  of  E and  k in  the  calculated  stability  boundaries. 

The  increase  in  amplitude  of  oscillation  with  increasing  alternating  voltage  was  large 
for  large  particles  — so  large  that  it  was  virtually  impossible  to  maintain  large  particles 
(>75  pm)  in  this  mode.  For  smaller  particles  the  increase  in  amplitude  was  less  marked;  particles 
of  20-pm  diameter  or  less  could  be  held  stably  at  any  point  in  the  instability  region.  These 
observations  aie  consistent  with  the  calculations  summarized  in  figure  10. 

Particles  of  18.5-pm  diameter  or  less  could  not  be  made  to  oscillate  horizontally  in 
the  double-period  mode.  For  these  particles,  k>  2.5 ; figure  7 shows  that  the  first  unstable 
region  does  not  extend  this  far  to  the  right,  so  double-period  vibration  is  not  to  be  expected. 
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When  the  alternating  voltage  was  just  below  that  required  to  produce  constant- 
amplitude,  double-period  oscillations,  and  the  laboratory  table  was  jarred,  a double-period 
oscillation  that  disappeared  in  a few  seconds  was  observed;  the  jarring  of  the  table  introduced 
a new  set  of  initial  conditions  and  new  values  of  A and  B in  equation  14,  so  that  the  amplitude 
of  the  fj  and  f2  functions  became  appreciable.  The  function  ^2  described  in  equation  18  no 
doubt  damped  out  almost  instanteously,  but  because  near  the  boundary  Ivl  = I 7 - K/2I  is  very 
small,  the  time  required  for  f | to  damp  out  and  become  negligible  was  long  enough  for  the 
oscillation  to  be  observed. 

One  calculated  result  not  observed  experimentally  is  shown  in  figure  12.  Calculation 
based  on  Oseen’s  Law  indicated  that  for  a 21-ftm  particle  (k  = 2.0  in  air  at  60  hertz),  as  the 
alternating  voltage  is  raised  above  the  value  required  to  produce  double-period  oscillation,  the 
amplitude  should  rise  monotonically  until  £ is  somewhat  above  the  upper  boundary  of  /?j  f/, 
and  then  drop  suddenly.  Qualitative  observations  of  20-ftm  particles  seemed  to  show  that  the 
amplitude  reaches  a peak  somewhere  in  the  Bi  [/  region  and  then  decreases  gradually  with  increas- 
ing voltage.  That  Oseen’s  modification  of  Stokes’  Law  is  itself  an  approximation  may  account  for 
this  discrepancy. 

B.  Experimental  Verifications. 

1 . Verification  of  the  Instability  Boundary  Curve. 

A 50%  solution  of  glycerine  and  water  was  placed  in  the  capillary  described  in  section  II. 
A single  particle  was  selected  and  inserted  into  the  chamber.  Figure  1 3 is  a photograph  of  a 
suspended  glycerine  particle. 


Figure  13.  A Suspended  Glycerine  Particle  80  /im  in  Diameter 
Focusing  of  the  backlight  by  the  particle  itself  causes  the  central  white  spot. 
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The  water  solution  was  used  because  it  accepted  more  charge  than  pure  glycerine.  The 
water  quickly  evaporated  leaving  an  almost  pure  glycerine  particle  that  was  highly  charged. 

Glycerine  was  chosen  as  the  calibrating  material  for  several  reasons.  First,  it  is  liquid 
at  operating  temperature.  As  a result,  the  particle  is  spherical  due  to  surface  tension.  Therefot^, 
the  diameter  d is  an  actual  value  rather  than  an  average  value,  making  the  calculation  of  k more 
accurate.  Second,  glycerine  is  soluble  in  water,  yet  not  highly  hygroscopic.  Hence  water,  used 
to  obtain  a more  highly  charged  particle,  could  be  more  easily  removed  by  placing  the  particle 
in  a low-humidity  chamber.  Finally,  the  vapor  pressure  of  glycerine  is  appropriate  for  the  time 
frame  required  to  make  these  measurements. 

The  ideal  situation  for  calibration  was  to  maintain  a decreasing  evaporation  rate  so 
that  the  rate  of  change  of  the  particle  diameter  remained  approximately  10  Mm/hr.  Glycerine 
particles  100  /rm  in  diameter  would  evaporate  far  too  slowly  at  room  temperature.  Therefore, 

T'  chamber  was  heated  at  the  beginning  of  the  run  and  then  slowly  cooled  during  the  time  the 
v;asurements  were  made.  Because  of  this  procedure  the  evaporation  rate  decreased  as  the 
temperature  decreased,  and  thereby  a fairly  constant  rate  of  change  of  the  particle  diameter  was 
maintained.  In  addition,  this  procedure  eliminated  thermaUy  induced  convection  currents. 

Typical  runs  started  with  a 1 00-/xm-diameter  particle  at  55®C  and  ended  with  a 
10-jum-diameter  particle  at  30°C.  Measurements  were  taken  at  various  time  intervals  during  the 
run.  The  procedure  consisted  of  several  steps.  First,  the  flow  of  dry  air  was  shut  off.  Next, 
the  diameter  of  the  particle  was  determined  using  the  calibrated  image-splitting  eyepiece.  The 
particle  was  then  balanced  by  varying  the  dc  voltage  until  the  two  images  of  the  particle,  formed 
by  operating  the  stroboscope  at  120  flashes  per  second,  coincided.  The  dc  voltage  was  then 
recorded  on  the  visicorder  and  oscilloscope.  Next,  the  ac  splitting  voltage  was  determined  and 
recorded  on  both  the  visicorder  and  the  oscilloscope.  Finally,  the  temperature  was  recorded  on 
the  visicorder. 

Because  all  of  the  previously  described  measurements  were  taken  within  1 minute, 
changes  in  the  particle  during  the  measurement  were  negligible. 

From  elementary  theory,  if  the  ac  voltage  is  increased  to,  or  becomes  greater  than, 
the  value  obtained  from  the  instability  boundary  curve  corresponding  to  a particular  particle 
diameter  and  density,  the  particle  should  begin  t'o  oscillate  at  one-half  the  driving  frequency  cu, 
with  increasing  amplitude.  The  amplitude  increases  until  the  particle  impacts  on  one  of  the 
electrodes.  As  a result,  it  appears  as  if  measuring  the  ac  splitting  voltage  requires  losing  the 
particle. 
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Because  of  nonlinear. effects  not  described  by  Stokes’  Law,  this  was  not  the  case.  At 
the  instability  boundary  the  particle  did  begin  to  oscillate  in  the  (j/2  mode.  However,  the 
amplitude  did  not  increase  with  time  until  the  ac  voltage  was  increased  beyond  the  stability 
boundary.  As  a result  the  ac  splitting  voltage  could  be  measured  without  losing  the  particle  by 
noting  the  sharp  transition  point  to  the  0)I2  mode. 


Figure  14  is  a graph  of  the  theoretically  calculated  instability  boundary.  This 
boundary,  as  shown  in  the  theoretical  section,  is  independent  of  chamber  design.  The  curve 
is  plotted  in  this  flgure  with  E and  d coordinate  axes.  Actually  the  boundary  is  a function  of 
E and  k as  previously  shown.  However,  the  diameter  d,  may  be  related  to  k by  the  equation 
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where  the  variables,  previously  described,  were  applied  to  the  case  of  a suspended  glycerine 
particle. 


Superimposed  *his  graph  are  experimental  data  points  computed  from  the  formula 


Qo^^ac 
^dc  ^^^dc 


using  the  actual  diameter  d in  micrometers. 

In  summary,  for  calculating  data  points  or  the  theoretical  curve,  g,  co,  rj,  and  p were 
known  constant  values.’  For  the  data  points,  and  d were  experimentally  determined, 

whereas  for  the  theoretical  curve  these  variables  were  allowed  to  assume  a continuous  range  of 
values. 


The  effect  of  chamber  design  was  included  by  means  of  the  parameters  by  Cac  and 
C(]q,  the  electric  field  linearity  constants.  Ideally,  both  parameters  could  be  calculated  theoreti- 
cally from  the  shape  of  the  electrodes.  In  fact,  Cac  was  easily  calculated  (as  shown  in  section  III.  A), 
but  the  calculation  of  failed  in  spite  of  many  approaches  attempted,  and  hence,  C^q  was 
experimentally  determined.  Section  III.  A shows  that  for  large  particles,  E limits  to  a constant 
value  of  0.464.  By  measuring  V^J^dc  wany  glycerine  particles  with  a diameter  greater  than 
1 1 0 pm  and  by  knowing  Cac>  ^ > and  w,  the  quantity  was  determined  to  be  17.7  ± 0.1 
meter  . 
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PARTICLE  DIAMETER 


FIELD  STRENGTH  PARAMETER  /: 


Figure  14.  Instability  Boundary  Curve 


2.  Rate  of  Instabilization  and  the  c>;/2  Boundary. 
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As  already  mentioned,  because  of  nonlinear  effects  not  described  by  Stokes’  Law,  a 
particle  maintains  stable  oscillation  beyond  the  instability  boundary.  After  the  boundary  is 
crossed,  the  particle  oscillates  with  increasing  amplitude  in  the  (j»l2  mode.  If  Stokes’  Law  were 
obeyed,  the  amplitude  would  continue  to  increase  and  the  boundary  would  be  a true  instability 
boundary.  However,  as  shown  in  the  theoretical  section,  nonlinear  drag  effects  become  impor- 
tant even  for  small  amplitudes  and  velocities  of  oscillation.  The  nonlinear  theory  predicts  that 
the  displacement  amplitude  increases  to  some  value  and  then  becomes  constant.  Therefore  the 
boundary  is  actually  a quasi-instability  boundary.  For  a given  percent  increase  in  E beyond  the 
boundary,  the  theory  also  predicts  that  larger  particles  have  larger  constant  amplitudes  than 
smaller  particles  in  the  same  state. 

Experimentally  it  was  found  that  for  larger  particles,  ac  voltage  increases  beyond  the 
splitting  voltage  cause  large-amplitude  changes.  This  is  not  true  for  smaller  particles.  Although 
quantitative  measurements  were  not  made,  this  phenomenon  is  completely  consistent  with  the 
theoretical  results  concerning  the  rate  of  instabilization. 

In  summary,  the  instability  boundary  is  actually  the  boundary  at  which  cj/2  frequency 
oscillation  begins  to  appear.  Because  of  nonlinear  drag  effects,  oscillation  at  this  frequency 
stabiUzes  at  an  amplitude  that  is  a function  of  the  value  of  E above  the  instability  boundary.  In 
this  sense  the  boundary  is  not  a true  instability  boundary. 

Although  the  cj/2  frequency  oscillation  begins  to  appear  at  this  boundary,  the  insta- 
bility boundary  is  not  the  point  at  which  the  aj/2  frequency  is  allowed.  As  shown  in  the 
theoretical  section,  the  co/2  boundary  is  an  entirely  different  curve.  However,  oscillation  in  the 
co/2  mode  is  damped  and  not  visible  until  the  instability  boundary  is  crossed. 

Experimentally  it  was  found  that  after  the  co/2  boundary  has  been  crossed,  but  before 
the  instabiUty  boundary  has  been  crossed,  the  particle  oscillates  only  in  the  co  mode  (if  unbalanced) 
On  the  other  hand,  if  the  particle  is  balanced  and  then  subjected  to  a vibration,  it  then  oscillates  in 
the  w/2  mode  until  the  motion  has  been  eUminated  by  damping.  Experimentally  it  was  also  found 
that  the  degree  of  co/2  mode  damping  lessens  as  one  approaches  the  instability  boundary;  this  is 
also  consistent  with  theoretical  results. 

3.  Passing  Through  the  First  Instability  Region. 

From  the  theoretical  description,  it  was  noted  that  the  first  instability  region  is  bounded 
and  separated  from  a second  instability  region.  When  the  value  of  E increases  from  zero,  a particle 
starts  in  a stable  state,  passes  through  an  unstable  region,  and  finally  ends  in  a stable  region  of 
motion. 
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It  appears  that  the  requirement  of  passing  through  the  primary  instability  region  would 
prohibit  observing  the  top  portion  of  the  first  instability  boundary.  However,  because  of  the 
quasi-instability  character  of  the  boundary  as  a result  of  nonlinear  drag  effects,  a particle  with  a 
diameter  of  20-/im  or  less  passes  through  this  region.  In  other  words,  the  constant  amplitude  of 
motion  in  the  cu/2  mode  did  not  exceed  the  distance  between  the  top  and  bottom  sheets  of  the 
hyperboloid  for  this  particle  diameter. 

When  smaller  particles  were  used,  the  top  portion  of  the  first  boundary  was  easily 
observed.  Starting  in  the  first  stable  region,  as  the  ac  voltage  was  increased,  the  particle  began 
to  oscillate  in  the  cj/2  mode,  indicating  that  the  first  boundary  had  been  reached.  As  higher 
voltages  were  applied,  larger  amplitude  oscillations  were  observed.  Finally  the  amplitude  decreased 
to  zero  indicating  that  the  top  portion  of  the  first  instability  region  had  been  crossed.  Further 
increase  of  the  ac  voltage  resulted  in  lower  values  of  the  restoring  force  as  the  second  instability 
region  was  approached.  At  the  second  boundary  there  was  no  restoring  force  at  all.  As  a result, 
this  boundary  was  never  experimentally  observed. 

1 

1 At  a particle  diameter  of  18.5  fxm  or  less  there  was  no  first  instability  boundary.  In 

this  case  it  was  observed  that  the  ac  voltage  could  be  varied  over  a large  range  of  values.  In  this 
I region  (k  > 2.4)  as  the  ac  voltage  was  increased  the  particle  would  move  toward  the  null  point 

iiV'  , until  a certain  voltage  was  attained.  If  a further  increase  of  voltage  were  applied,  the  particle 

would  then  move  away  from  the  null  point.  This  indicated  that  the  net  restoring  force  reached 
( a maximum  as  the  ac  voltage  was  increased  with  this  maximum  being  dependent  on  the  particle 

size  as  shown  in  figure  10. 

4.  Phase  Angle. 

As  mentioned  in  the  introduction  and  described  in  detail  in  section  III. A,  the 
phase  angle  between  the  ac  voltage  oscillation  and  the  particle  oscillation  was  the  basic  parameter 
in  accounting  for  the  nonzero  time-averaged  force  toward  a point  in  space  experienced  by  the 
particle.  The  value  of  the  phase  angle  varied  with  different-sized  particles  and  could  have  been 
used  to  measure  particle  diameter. 

Experimental  measurement  of  the  phase  angle  was  accomplished  by  flashing  the 
stroboscope  at  1 20  hertz.  The  time-delay  attachment  to  the  stroboscope  was  varied  until  the 
two  images  of  the  particle  coincided,  at  which  point  the  stroboscope  was  flashing  in  phase  with 
the  particle  motion.  The  phase  angle  was  then  determined  from  the  displaying  of  the  ac  voltage 
and  the  stroboscope  flash  on  a two-channel  oscilloscope. 

As  shown  in  section  111. A (figure  11),  actual  phase  angle  measurements  made  for 
the  bihemispheric  system  with  a 50%  solution  of  water  and  ethylene  glycol  were  in  excellent 
agreement  with  the  calculated  values.  However  similar  phase  angle  measurements  were  not  made 
for  the  bihyperboloidal  system  described  in  this  report. 


5.  Particle  Discharge. 


Particle  discharge  has  been  investigated  and  reported  by  others*.  In  theory  particle 
discharge  should  occur  when  the  electrostatic  force,  which  increases  during  evaporation  of  the 
droplet  as  the  charge-to-mass  ratio  increases,  equals  the  surface  ten.sion  force.  This  condition  is 
commonly  known  as  the  Rayleigh  stability  limit.  The  nature  of  particle  discharge  has  been 
described  as  the  ejection  of  several  smaU  highly-charged  droplets  comprising  about  30%  of  the 
electric  charge  and  5%  of  the  mass  of  the  original  droplet. 

In  practice  particle  discharge  was  noted  by  the  sudden  disappearance  of  the  droplet. 
If  the  ac  voltage  was  high  enough  to  contain  the  particle,  the  particle  could  be  returned  to  the 
field  of  view  by  increasing  the  dc  voltage.  The  illumination  from  the  stroboscope  was  not 
sufficient  to  distinguish  any  small  ejected  droplets;  however,  the  electrical  measurements  on  the 
large  droplet  after  discharge  indicated  a 20%  to  40%  loss  of  charge.  The  5%  decrease  in  mass 
could  not  be  measured  accurately  because  of  the  time  required  for  rebalancing  the  particle  to 
compensate  for  the  charge  loss. 


In  the  chamber  at  the  time  of  discharge,  the  forces  should  be  related  by 

_ 3.768pV 

Qc^^'dc^ 


where  a is  the  surface  tension  and  the  right-hand  side  is  the  electrostatic  force.  Calculations  for 
several  droplets  were  made  from  the  data  just  prior  to  discharge  of  the  droplet.  These  results 
varied  considerably  from  the  measured  surface  tension  of  the  bulk  liquid  (63  dynes/cm)  and 

ranged  in  value  from  20  to  SO  dynes/cm.  No  apparent  reason  for  this  wide  range  of  values  has 
been  determined. 

6.  Charge  Decay  Due  to  Cosmic  Ray  Ionization. 

One  source  of  interest  and  possible  error  was  the  rate  of  loss  of  charge  due  to  cosmic 
ray  ionization.  A polyethylene  bead,  66.7  pm  in  diameter,  was  placed  in  the  chamber  and  sus- 
pended for  68  hours.  Loss  of  charge  was  indicated  by  an  increase  in  the  dc  balancing  voltage. 
The  measured  rate  of  increase  was  approximately  1 .76  X I O"^  volts/second. 

With  dry  air  circulating  through  the  chamber,  the  atmosphere  was  free  of  contaminants. 
A logical  explanation  of  the  discharge  seemed  to  be  ionization  caused  by  cosmic  ray  radiation. 


The  rate  of  ion-pair  production  at  sea  level  is  1.5  to  2 ion  pairs  per  sec-cm^.‘®  These 
ions  are  small  and  have  a mobility  of  about  2 cm^/(v-«ec).  Thus  they  can  only  exist  within  the 
chamber,  impinging  on  the  chamber  walls  or  on  the  particle,  for  less  than  one^lf  of  a cycle  of 
the  ac  voltage. 

Calculations  showed  that  the  charge  on  the  polyethylene  particle  was  approximately 
1.8  X 10~12  coulombs.  It  was  found  that  the  electric  field  due  to  the  particle  equaled  the  ac  field 
at  a distance  of  O.S  cm  from  the  particle.  Based  on  this  fact  and  the  design  of  the  chamber,  one 
would  expect  that,  during  one-half  of  a cycle,  ions  contained  in  a disk  2.50  cm  in  diameter  and 
1 cm  in  length  and  ions  contained  in  a cylinder  3.50  cm  in  length  and  1 cm  in  diameter,  both 
coaxial  with  the  symmetry  axis  of  the  chamber,  would  impinge  on  the  particle.  Thus  the  rate  of 
discharge  would  be  approximately  2.9  to  3.9  electronic  charges  per  second.  (Note  that  one-half 
of  the  ion  pairs  produce  a discharge  for  one-half  of  a cycle.) 


From  the  relation 


^ = jni.  Fd  2 
dt  Cdc 


dt 


where  all  the  quantities  have  been  previously  described,  the  rate  of  increase  of  f^dc  found  to 
correspond  to  a discharge  rate  of  3.8  electronic  charges  per  second. 
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APPENDIX  A 


DERIVATION  OF  COEFFICIENTS 


The  relation  between  equation  1 5 


and  the  Mathieu  equation 


f"(r)  + Kt'  + (^i  + £■  cos  r)i’(r)  = 0 


(t)  + (a  -2q  cos  2t  )J:(t)  = 0 


is  described,  for  instance,  in  reference  1 1 . 


T = 2t 


f = e-KT  j 


f'(T)  = ^ e->^^  (f'  -<^T) 

f"(T)  = j-  [p'  - 2/ct'  + ic^f] 
4 


Substituting  equation  A-2  into  equation  1 5 yields 


+ (4£'j  - + 4^  cos  2T)f(T)  = 0 


This  is  equation  A-1  where 


= 4£j  - ic2 


q = -2£ 


as  described  in  reference  1 1 . 


jgg[SQUIQt  PAflB  MOT  HU®) 


The  solutions  to  equation  A-1  are  of  the  form 


n=-oo 

oo 

«=-0O 

where  y takes  one  of  the  forms 

y = a 
y = a + i 


(A-t) 


(A-5) 


y = 0i 

where  a and  0 are  real  and  a is  nonnegative. 

From  equations  17,  A-4,  and  A-5  it  is  evident  that  the  solutions  of  equation  15  are 


n=  -«» 

OQ 

n=-°° 


where  y may  be  of  the  form  y = a,  y = a + or  y = 0i  for  real  a and  0. 


The  values  of  a and  0 as  functions  of  a and  q have  been  calculated  and  are  displayed, 
e.g.,  in  reference  1 1,  as  Ince’s  chart.  It  is  possible  to  take  values  from  Ince’s  chart  and  use  them 
for  the  chart  shown  in  figure  7 via  equations  1 7 and  A-3.  In  this  investigation,  however,  calcu- 
lations were  made  directly  from  equations  1 5 and  1 8 by  the  method  of  continued  fractions 
described  in  references  1 1 and  1 2. 


The  solution  to  the  inhomogeneous  form  of  equation  1 5 given  by  equation  1 3 

f"  + Kf'  + (fj  + E cos  r)f  + h = 0 (13) 
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Then  equation  A-7  can  be  written  as 


whence 


If  r„  is 


iE^-n^+inK)■^  j-±  + |r„+i  =0 


(A-10) 


j (n2-£i-/n(c)-r„+i 


(A-11) 


fn  = 


(A-12) 


equation  A-10  yields 


= 


J-  (n^-Ei  +inKh  r„^l 


(A-13) 


It  is  evident  that  for  any  value  of  «q,  r^Q  can  be  expressed  as  a continued  fraction  in 
terms  of  where  n > kq. 

To  assume  that  the  series  representing  p(t)  converges,  we  require  that  for  some  value  of 
hq,  \r„  I and  l7„  I be  less  than  some  upper  bound  if  « > mq.  Then  for  large  enough  mj 


+ 1 1 « 4 /hkI 


and 


'•«!  + ll«  +/„^| 


(A-14) 


For  n = nj,  equation  A-1 1 may  be  written 

1 


f-  - i/ijK) 


^«1  = 


E 


(A-1 5) 
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Using  these  values  of  r„j  and  ^nj,  may  be  calculated  from  equa- 

tion A-1 1.  In  the  same  way  values  of  and  7„  may  be  calculated  for  successively  lower  values 
of  n,  until  rj  and  7j  have  been  found.  Equations  A-1 5 show  that,  since  E,n.  Ei,  and  7 are 
real,  *‘1'****°***  ^ show  that  if  r;,+  i * 7„*+i,  then  • TJi*.  Hence 


= 'n* 


(A-1 6) 


for  all  /I,  and  only  the  r„  need  be  calculated. 

From  equations  A-9,  A-1 2,  and  A-1 6 


aj  = rjao 


fl-1  = 


whence  equation  A-8  can  be  written 


tfQffj  +£>?(ri)l  + A = 0 


(A-17) 


‘*0'"  £1  +£i?(ri) 

where  Rir^)  is  the  real  part  of  rj.  The  value  of  flg.  th®  ™®an  position  of  the  particle  with  respect 
to  the  null  point,  is  real,  as  it  must  be,  being  a physical  quantity.  Then  /i(ai)  = and 

equation  A-17  can  be  written 

ERiai)  + h + oqEi  = 0 

This  equation  is  essentially  equation  21  reached  in  a different  way. 

Starting  from  equation  A-1 8,  successive  coefficients  of  p(t)  may  be  calculated  as 


3n+l 
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In  actual  calculation,  the  starting  value  of  n need  not  be  so  high  that  equation  A- 14 
applies;  the  error  in  decreases  very  rapidly  as  n decreases.  For  £ < 10,  the  value  of  need 
exceed  by  no  more  than  10  the  highest  value  of  n for  which  an  accurate  value  of  a„  is  required. 

In  order  to  determine  the  value  of  y in  equation  18,  the  expression  for  fi  is  substituted 
into  equation  IS  yielding 

oo 

+ £j  + 2lny)b„  + + 6„.l  )I  * 0 


whence 


+ b„.i)  (n^  + ^K^-y^-E\  - 2iny)b„  = 0 (A-19) 


Let  r„  and  be 


" bn.l 


respectively. 


Then  equation  A-19  is 


- +'-n+l  = f - y^-Ei-2iny) 


-|(«2  + 1x2-  72.  £j . 2my)- r„+. 


” -y^-E\  + 2iny)  - r„ 


(A-20) 


(A-21) 


As  before,  r„+i  is  regarded  as  negligible  for  some  sufficiently  large  value  of  n,  and  r„ 
7ff  are  calculated  for  successively  lower  n.  For  n = 0,  equation  A-19  is 


(A-22) 


1 

I 

J 
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In  general,  equation  A'22  is  not  true.  If  any  three  of  the  four  quantities  E,  E^,  k,  and 
7 are  known,  a value  of  the  fourth  quantity  must  be  sought  for  which  equation  A-22  is  satisfied. 
In  this  investigation  the  entities  to  be  determined  were  boundaries  between  various  regions  and 
curves  of  constant  value  of  p,  where 


p = a ~~K 
2 


in  regions  where  7 # 0. 


On  these  curves 


y ~ p + ~K 
2 


7 = *'  + + 4-/ 


On  periodicity  boundaries  between  regions  where  7 = 0 and  7 = /Ji,  it  is  evident  that 
unless  7 is  a discontinuous  function  of  k and  E, 

7 = 0 

On  periodicity  boundaries  between  7 = 0 + -j/  and  7 = /Ji 

On  stability  boundaries  a = 

7 = |»c 

or 

7 = 4-' 

2 2 

Except!  for  the  calculations  leading  to  figure  8,  = 0.  In  any  case  £]  was  taken  as 

a known  quantity.  Then,  depending  on  the  part  of  the  curve  being  calculated,  either  k was  fixed 
and  corresponding  values  of  E were  sought,  or  vice  versa. 


On  periodicity  boundaries  where  7 = 0,  equations  A-20,  A-21,  and  A*22  are 
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(A-23) 


” ' 'n+l 


It  i>  evident  that  and  are  real  and 


so  that  equation  A-23  is 


'■1  = 


On  periodicity  boundaries  where  7 = -ji 


I ((n+^)^  + -'■«+! 


|r[(«“)2  + - ’n+l 


“ '’n-l 


and  equation  A-23  is 


On  stability  boundaries  where  y = -jK, 


" |r(n2-£j-  lnic)-r„+, 
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SO 


1 


•A* 


-J  T 

.£.  +m»c)  -r„+i 
E 


fr.  = »■« 


and 


_^1 


Rir,)=-^ 


On  stability  boundaries  where  7 = -^k  + 


= 


^[in+^9- - E^- iin-¥^)K\- r„+l 


1 


|r[(«-^)2-£i  '•«+! 


On  lines  of  constant  v where  y = v + -^k, 


'■«  = 


1 


-2-  [«2  . i;(i;  + »c)-  £1  - l>l(2»'  + K)]  - r„4.] 
£ 


r„  = 


-2-  [rt^  - t'fi'  + k)  - £1  + + k)1  * ^n+1 

£ 
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Rir^)  = + k)  + £j  1 


On  lines  of  constant  v where  y = v + -i-ic  + -i-i, 

2 2 


” -|:f  (n+-y)^-  «'(»'  + k)-  - i(n+-^)(2v  + «))-  r„+i 


«'(«'  + »<)-£i  +i(n  --1)(2i;  + k)1  -r„+i 


-F  - r* 

- '^n-l 


(A-24) 


ri  ''O  ^ £■  t "4  + >')  - 1 


(A-25) 


For  each  of  these  cases  the  scheme  of  calculation  v/as  as  follows: 

1.  The  value  of  fj  was  fixed,  usually  at  =0. 

2.  For  the  case  of  lines  of  constant  v,  the  value  of  i>  was  fixed. 

3.  The  value  of  k (or  in  some  cases  the  value  of  £)  was  fixed. 

4.  Using  a guessed  value  of  £ (or  «),  r j and,  if  necessary,  tq  were  computed  by  repeated 
application  of  equation  A-20,  starting  from  some  sufficiently  high  value  of  n.  Note  that  equa- 
tion A-21  is  not  us  d in  the  computation,  being  required  only  to  establish  the  relation 

5.  The  difference  between  the  two  sides  of  equation  A-25  was  determined. 

6.  A new  guess  at  £ (or  k)  was  made,  usually  by  interpolation. 

7.  When  a value  of  £ (or  k)  was  found  for  which  equation  A-25  was  sufficiently 
nearly  satisfied,  a new  k (or  £)  was  fixed  and  steps  4 through  6 repeated. 

The  curves  shown  in  figure  1 2 were  calculated  by  writing  the  Taylor  expansion  for 
i and  i'  from  equation  1 1 . 
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f"  = -(«  + «!«)?'  - {£  cos  T + £j  )f  - /« 

f'"  = -(k  + Kjw)f" . («!«’  +£  cos  r + )f'  + sin  r 

{•(t  +Af  ) = f(r)  + Arf'(r)  + 4" 

2 o 

r (T  + A?  ) = f'(r)  + A?  f "(r)  + ^ A?  2 f "'(r) 

where  At  is  the  amount  by  which  r is  incremented,  usually  7r/12S  (that  is,  cycle  of  cos  t), 
and  u is  the  magnitude  of  the  velocity  as  represented  by  I f ' L In  these  computations  u was 
taken  as  If'  I and  the  contribution  of  the  horizontal  component  of  the  velocity  was  ignored, 
even  though  where  the  horizontal  unstable  region  overlaps  the  vertical  stable  region,  the  hori- 
zontal contribution  is  not  negligible.  Most  of  the  region  in  k-E  space  where  there  is  instability 
in  both  directions  has  been  hitherto  inaccessible  to  experimental  observation  in  any  case.  Near 
the  nose  of  the  first  instability  curve  where  observation  in  the  first  unstable  region  is  easy,  the 
horizontal  contribution  damps  out. 

Computations  were  carried  out  on  Mathetron  calculators  and  on  a Univac  1 1 08  computer. 
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LIST  OF  SYMBOLS 


integration  constant  in  solution  to  dimensionless  equation  of  motion 

coefficients  in  expression  for  potential  d>  in  cylindrical  coordinates 

coefflcients  in  expression  for  potential  $ in  spherical  coordinates 

coefficient  in  Mathieu’s  equation 

coefficients  of  Fourier  expansion  of  p(t)  in  terms  of 

integration  constant  in  solution  to  dimensionless  equation  of  motion 

ratio  of  gradient  of  vertical  component  of  electric  field  at  null  point  to  potential 
difference  between  electrodes 

constant  given  by  Mg  | 

constants  used  to  define  electrode  hyperboloids  | 

ratio  of  field  strength  at  center  of  chamber  to  potential  difference  between  sheets 
of  hyperboloid 

coefficients  of  Fourier  expansion  of  p(t)  in  terms  of  cos  {nr  ± <l>„) 
diameter  of  particle 

peak  value  of  field  strength  due  to  alternating  potential  difference  between 

electrodes  i 

I 

field  strength  due  to  constant  potential  difference  between  electrodes  J 

Ej  + E cos  ut 

field  strength  due  to  potential  difference  between  sheets  of  two-sheet  electrodes 


e 2.718281828... 

E coefficient  of  f cos  w/  in  dimensionless  equation  of  motion  (due  to  E) 

£j  coefficient  of  f in  dimensionless  equation  of  motion  (due  to  Ej) 

g acceleration  due  to  gravity 

g magnitude  of  g 

h constant  term  in  dimensionless  equation  of  motion  (due  to  gravity  and  Ejjj.) 


1 


ll 

I., 


K 


D 


rn 


P 

<i 

Q 

R 

R(X) 

Re 

^0 

R„S 

R„U 


’’0 

''"o 


u 

V 


dc 


first  region  where  y = fii 
imaginary  part  of 

coefficient  describing  aerodynamic  drag 

unit  vector  along  vertical  axis  (r  axis) 

mass  of  particle 

summation  index 

nth  Legendre  polynomial 

periodic  factor  in  solution  to  equation  of  motion 

periodic  solution  to  inhomogeneous  equation 

charge  on  particle 

coefficient  in  Mathieu’s  equation 

position  of  particle 

real  part  of  X 

real  part  of 

Reynolds  number 

lowest  7 = 0[  region 


stable  part  of  7 = a or  7 = a region 
unstable  part  of  7 = a or  7 = a + region 
radial  coordinate  in  cylindrical  coordinate  system 
radial  coordinate  in  spherical  coordinate  system 
inner  radius  of  one-sheet  hyperboloidal  electrode 
unit  radial  vector  in  cylindrical  coordinate  system 
time 

magnitude  of  particle  velocity  divided  by  particle  diameter 

peak  value  of  alternating  part  of  potential  difference  between  electrodes 

constant  part  of  potential  difference  between  electrodes 

potential  difference  between  sheets  of  two-sheet  electrode 
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y 

z 

^0 

a 

a 

0 

? 

7 

7 

fi’  h 

F],  f2 

V 

0 

K 

V 


horizontal  cartesian  coordinates 
vertical  cartesian  or  cylindrical  coordinate 

half-height  of  hyperboloidal  chamber  (distance  from  null  point  to  nose  of  two- 
sheet  hyperboloidal  electrode) 

real  part  of  7 for  equation  of  motion 

real  part  of  7 for  Mathieu’s  equation 

imaginary  part  of  7 for  equation  of  motion 

imaginary  part  of  7 for  Mathieu’s  equation 

exponential  coefficient  in  solution  to  equation  of  motion 

exponential  coefficient  in  solution  to  Mathieu’s  equation 

solution  to  dimensionless  equation  of  motion 

two  independent  solutions  to  dimensionless  equation  of  motion  for  vertical  motion 
solutions  to  Mathieu’s  equation 
viscosity  of  air 

equatorial  angle  in  cylindrical  coordinate  system 

polar  angle  in  spherical  coordinate  system 

drag  term  coefficient  in  dimensionless  equation  of  motion 

velocity  squared  drag  term  coefficient  in  dimensionless  equation  of  motion 


solutions  to  dimensionless  equations  of  motion  for  horizontal  motion 
3.14159265  . . . 


p density  of  air 

Pj  density  of  particle  material 

T (independent  variable  in  dimensionless  equation  of  motion) 

? independent  variable  in  Mathieu’s  equation 

TQ  value  of  T at  coincidence  point 

4>  potential  function 
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■I 


a 

OJ 

asterisks 

dots 

primes 


lag  angle  for  nth  term  of  series  representing  p(r)  in  terms  of  cosines 
n<l>i  - 

surface  tension  of  particle  material 

angular  frequency  of  alternating  potential  difference  between  electrodes 
complex  conjugate 

differentiation  with  respect  to  time  (f) 
differentiation  with  respect  to  r 
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RESEARCH  DIVISION 

Atm:  DRDAR<LB  I 

Attn:  DRDAR-CLB-B  1 

Atm:  DRDAR-CLB-C  1 

Atm:  DRDARCLB-P  1 

Atm:  DRDAR-CLB-R  1 

Atm.  DRDARCLB-T  1 

Atm:  DRDAR-CLB-TE  1 

SYSTEMS  ASSESSMENTS  OFFICE 

Atm:  DRDAR-CLY-A  1 

Atm:  DRDAR-CLY-R  1 


DEPARTMENT  OF  DEFENSE 
Administrator 

Defense  Documentation  Center 

Atm:  Accessions  Division  (DDC-TC)  12 

Cameron  Station 
Alexandria,  VA  22314 

Director 

DefenK  Intelligence  Agency 

Atm:  DB-4G1  1 

Washington,  DC  20301 


US  ARMY  MATERIEL  DEVELOPMENT  AND 
READINESS  COMMAND 

Commander 

US  Army  Materiel  Development  and  Readiness  Command 
Atm:  DRCLDC  1 

Atm:  DRCSF-P  1 

5001  Eisenhower  Ave 
Alexandria,  VA  22333 

Office  of  the  Project  Manager  for  Chemical  Demilitarization 
and  Installation  Restoration 

Atm:  DRCPM-DR-T  2 

Aberdeen  Proving  Ground,  MD  21010 

Human  Engineering  Laboratory  HFE  Detachment 

Atm;  DRXHE-EA  1 

BuOding  E3220 
APG-Edgewood  Area 

Commander 

US  Army  Foreign  Science  4 Technology  Center 

Atm:  DRXST-CX2  I 

220  Seventh  St.,  NE 
Cliarlotteiville,  VA  22901 
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Commander 

US  Army  Missile  Research  and  Development  Command 
Redstone  Scientific  Information  Center 
Attn:  DRDMI-TBD 
Redstone  Arsenal,  AL  3S809 


Commander 

Rocky  Mountain  Arsenal 
Attn:  SARRM-QA 
Attn;  SARRM-MD 
Commerce  City,  CO  80022 


Director 

DARCOM  Field  Safety  Activity 
Attn.  DRXOS-C 
Charlestown,  IN  471 1 1 


' Commander 
Pine  Bluff  Arsenal 
1 Attn;  SARPB-ETA 

Pine  Bluff,  AR  71611 


Commander 

US  Army  Materiel  Development  and  Readiness  Command 
Installations  and  Services  Activity 
Attn:  DRCIS-RI-IU 
Rock  Island,  IL  61299 

US  ARMY  ARMAMENT  RESEARCH  AND 
DEVELOPMENT  COMMAND 

. mmander 

US  Army  Armament  Research  and  Development  Command 


Attn:  DRDAR-LCE-M  1 

Attn:  ERDAR-LCH  1 

Attn:  DRDAR-LCU  1 

Attn:  DRDAR-SER  1 

Attn:  DRDAR-TSS  2 

Dover,  NJ  07801 

Director 

Ballistic  Research  Laboratory 

Attn:  DRDAR-TSB-S  1 

Building  328 


Aberdeen  Proving  Ground,  MD  21005 

CDR,  APG 
USA  ARRADCOM 

Attn:  DRDAR-GCL  1 

Aberdeen  Proving  Ground,  MD  21010 

US  ARMY  ARMAMENT  MATERIEL  READINESS  COMMAND 


US  ARMY  TRAINING  & DOCTRINE  COMMAND 
Commandant 

] US  Army  Infantry  School 
Attn:  NBC  Division 
Fort  Benning,  GA  31905 

Commandant 

US  Army  Missile  & Munitions  Center  A School 
Atm:  ATSK-CD-MD 
Atm:  ATSK-DT-MU-EOD 
Redstone  Arsenal,  AL  35809 

Commandant 

US  Army  Military  Police  School/Training  Center 
Atm:  ATZN-CDM 
Atm:  ATZN-TDP-C 
Fort  McOellan,  AL  36205 

Commander 

US  Army  Infantry  Center 
Atm:  ATSH-CD-MSC 
Fort  Benning,  GA  31905 

Conun  andant 

US  Army  Ordnance  A Chemical  Center  A School 
Atm:  ATSL-CL-CD 
Aberdeen  Proving  Ground,  MD  21005 

US  ARMY  TEST  A EVALUATION  COMMAND 


Commander 

US  Army  Armament  Materiel  Readiness  Command 
Atm:  DRSAR-ASN 
Atm:  DRSAR-IMBC 
Atm:  DRSAR-SF 
Rock  Island,  IL  61299 

CDR,  APG 
USA  ARRCOM 
Atm:  SARTE 

Aberdeen  Proving  Ground,  MD  21010 
Commander 

US  Army  Dugway  lYovlng  Ground 
Atm:  Technical  Library,  Docu  Sect 
Dugway,  UT  84022 


Commander 

US  Army  Test  A Evaluation  Command 
1 Atm:  DRSTE-FA 

1 Aberdeen  Provmg  Ground,  MD  21005 

1 

Commander 

US  Army  Cold  Regions  Test  Center 
Atm:  STECR-TD 
APO  Seattle,  WA  98733 
I 

DEPARTMENT  OF  THE  NAVY 

Chief  of  Naval  Research 
Atm:  Code  443 
1 800  N.  Quincy  Street 

Arlington,  VA  22217 


in 


Copies 
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1 
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1 


1 


1 
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Commander 

HO  AFISC/SEV 

1 

Naval  Explosive  Ordnance  Disposal  Facility 

Norton  AFB,  CA  92409 

Attn:  Army  Chemical  Officer,  Code  604 

1 

Indianhead,  MD  20640 

Commander 

NORAD  Combat  Operations  Center/DBN 

Cheyenne  Mtn  Complex,  CO  80914 

1 

Naval  Surface  Weapons  Center 

White  Oak  Laboratory 

OUTSIDE  AGENCIES 

Attn:  Tech  Lib  & Info  Svcs  Br 

1 Battelle,  Columbus  Laboratories 

Silver  Spring,  MD  20910 

Atm:  TACTEC 

505  King  Avenue 

1 

Commander 

Nuclear  Weapons  Training  Group,  Atlantic 

Columbus.  OH  43201 

Naval  Air  Station 

Director  of  Toxicology 

1 

Attn:  Code  21 

1 National  Research  Council 

Norfolk,  VA  23511 

2101  Constitution  Ave,  NW 

Washington,  DC  20418 

Chief,  Bureau  of  Medicine  & Surgery 

1 

Department  of  the  Navy 

Washington,  DC  20372 

ADDITIONAL  ADDRESSEES 

US  Public  Health  Service 

Commander 

Room  17A-46  (CPT  Osheroff) 

1 

Naval  Weapons  Center 

5600  Fishers  Lane 

Attn;  A.  B.  Calloway/Code  3171 

1 Rockville,  MD  20857 

Attn;  Technical  Library/Code  233 

1 

China  Lake,  CA  93555 

Commander 

US  Army  Environmental  Hygiene  Agency 

US  MARINE  CORPS 

Director,  bevelopment  Center 

Atm:  Librarian,  Bldg  2100 

Aberdeen  Proving  Ground,  MD  21010 

1 

Marine  Corps  Development  & Education  Command 

Commander 

Attn:  Fire  Power  Division 

1 DARCOM,  STITEUR 

Quantico,  VA  22134 

Atm.  DRXST-STl 

Box  48,  APO  New  York  09710 

1 

DEPARTMENT  OF  THE  AIR  FORCE 

Commander 

HO  Foreign  Technology  Division  (AFSC) 

US  Army  Science  & Technology  Center-Far  East  Office 

1 

Attn:  PDRR 

Wright-Patterson  AFB,  OH  45433 

1 APO  San  Francisco  96328 

HQDA  DASG-RDZ  (SGRD-PL) 

1 

Commander 

Aeronautical  Systems  Division 

WASH  DC  20314 

Atm;  ASD/AELD 

I Commander 

Wright-Patterson  AFB,  OH  45433 

USEUCOM 

Attn:  EOS-O/LTC  Jainei  H.  Alley 

1 

AMRL/MEB 

Wright-Patterson  AFB,  OH  45433 

1 APO  New  York  09128 

HO.  USAF/SGPR 

Forrestal  Bldg 

WASH  DC  20314 

1 

HO  USAF/RDPN 

WASH  DC  20330 

1 
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